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PREFACE 

In  the  past  several  years,  the  National  Transportation  Safety  Board  has 
identified  wind  shear  as  a factor  or  probable  cause  in  several  aircraft  acci- 
dents. The  most  significant  of  the  wind  shear-related  accidents  have  been 
Iberia  Flight  933  at  Boston  (a  DC- 10),  Eastern  Airlines  Flight  66  at  New  York 
(B-727),  and  Continental  Airlines  Flight  426  at  Denver  (B-727).  Wind  shear  is 
a gradient  of  wind  velocity  with  respect  to  altitude  above  the  surface.  The 
accidents  have  prompted  a search  for  techniques  and/or  instrumentation  for 
the  remote  measurement  of  wind  as  a function  of  altitude  for  the  purpose  of 
identifying  regions  of  wind  shear  and  determining  the  magnitudes  and  direction 
of  the  shear.  A primary  candidate  for  the  remote  measurement  of  wind  is  the 
laser  Doppler  velocimeter.  The  test  described  herein  was  designed  to  validate 
the  laser  Doppler  velocimeter  as  an  instrument  for  the  remote  measurement 
of  wind  and  to  evaluate  alternative  data  processing  algorithms.  This  report 
describes  the  principles  of  wind  measurement  by  laser  Doppler  velocimetry 
and  presents  the  results  of  the  validation  tests. 
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1.  INTRODUCTION 

1.1  BACKGROUND 

Significant  effort  is  currently  being  devoted  to  the  development  of  in- 
strumentation to  sense  atmospheric  flow  phenomena  remotely.  Some  of  the 
avenues  being  pursued  are  active  and  passive  acoustic  sensors,  optical 
sensors,  and  radio  methods.  A useful  survey  of  such  methods  is  presented 
in  Ref.  1.  Two  advantages  of  remote  sensors  are  that  flow  conditions  can 
be  ascertained  in  regions  of  space  where  it  would  not  be  convenient  to  locate 
instrumentation  hardware,  and  no  interference  with  the  flow  at  the  point  of 
interest  is  introduced  by  their  use.  The  Laser  Doppler  Velocimeter  (LDV) 
is  a particularly  attractive  device  for  remote  sensing  of  atmospheric  phe- 
nomena. In  the  LDV  system,  the  laser  radiation  backscatter ed  by  moving 
particulates  in  the  atmosphere  is  used  to  determine  the  velocity  of  the  flow. 
Since  it  is  possible  to  direct  the  laser  focal  volume  at  a selected  sequence 
of  points  in  space,  data  from  a scanning  LDV  system  can  be  used  to  deter- 
mine the  velocity  field  rapidly  and  over  a range  of  altitudes.  A CO^  laser 
Doppler  velocimeter  system  has  the  following  advantages  over  other  remote 
sensing  techniques:  (1)  the  sensing  volume  can  be  varied  with  ease  as  only 
optic  pointing  and  focusing  operations  are  involved;  (2)  the  ambient  aerosol 
provides  a sufficient  scattering  target;  and  (3)  the  sensing  mechanism  is  non- 
mechanical which  results  in  the  potential  for  a high-frequency  turbulence 
sensor. 

The  feasibility  of  utilizing  an  LDV  system  for  the  remote  sensing  of 
low  altitude  winds  and  for  the  detection  and  tracking  of  aircraft  wake  vortices 
has  been  demonstrated  (Refs.  2 through  5).  However,  the  development  of  an 
effective  LDV  system  for  monitoring  wind  and  wind  shear  required  a further 
refinement  and  application  of  the  technology. 
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1.2  PROGRAM  OBJECTIVES 

The  primary  objectives  of  the  Table  Mountain  test  were  to  establish 
the  validity  of  the  Velocity  Azimuth  Display  (VAD)  for  three-dimensional  wind 
measurements  for  a range  of  wind  speeds  and  to  assess  the  validity  of  se- 
lected VAD  data-processing  techniques.  Another  objective  was  to  investigate 
the  effect  of  averaging  time  on  the  accuracy  of  VAD  measurements. 

1.3  REPORT  FORMAT 

A brief  description  of  the  mobile  laser  Doppler  system  is  presented 
in  the  following  sections.  A more  complete  description  is  presented  in 
Ref.  5.  A description  of  the  LDV  is  presented  in  Section  2 followed  by  a 
description  of  the  computer  software  algorithms  in  Section  3.  The  field 
tests  at  the  Table  Mountain  site  are  described  in  Section  4 along  with  analy- 
sis of  the  VAD  data  and  comparisons  with  the  meteorological  tower  data. 
Section  5 presents  the  conclusions  and  recommendations. 
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2.  LASER  DOPPLER  SYSTEM  DEVELOPMENT 
2.1  SYSTEM  DESCRIPTION 

An  LDV  remote  wind  sensor  senses  air  movement  by  measurement  of 
the  Doppler  frequency  shift  of  laser  radiation  backscattered  by  the  atmos- 
pheric aerosol.  An  instrument  must  incorporate  means  to  transmit  the  laser 
radiation  to  the  region  of  interest,  collect  the  radiation  scattered  from  the 
atmospheric  aerosol,  and  to  photomix  on  a photodetector  the  scattered  radia- 
tion and  a portion  of  the  transmitted  beam.  The  d tference  between  the  trans- 
mitted frequency  and  the  returned  frequency  is  the  Doppler  shift  frequency. 

The  Doppler  frequency  shift  signal  is  generated  at  the  photodetector  and  is 
directly  proportional  to  the  magnitude  of  the  wind-velocity  component  in  the 
direction  of  the  line-of-sight  of  the  laser  beam.  This  velocity  component  is 
hereafter  called  the  line-of-sight  velocity.  The  magnitude  of  the  Doppler  shift, 
Af,  is  given  by 

Af  = — I v| cose, 

|A 

where 

V = the  velocity  vector)  in  the  region  being  sensed, 

9 = the  angle  subtended  by  the  velocity  vector  and 
the  optic  system  line-of-sight,  and 

X = the  laser  radiation  wavelength  (10.6  microns 
for  CO^  laser). 

A Doppler  shift  of  188  kHz  results  per  m/sec  of  line-of-sight  velocity  com- 
ponent. Measurement  of  the  Doppler  shift  frequency,  Af,  yields  directly  the  line- 
of-sight  velocity  component,!  Vl  cos0.  Some  typical  advantages  of  the  laser 
Doppler  method  are;  (1)  the  Doppler  shift  is  a direct  absolute  measure  of  the 
line-of-sight  velocity  (for  example,  a hot  wire  yields  wind  speed  via  a cooling 
effect  on  the  wire);  (2)  the  ease  with  which  the  sensing  volume  can  be  varied 
(only  optics  pointing  and  focusing  operations  being  involved);  (3)  the  ambient 

-Jf 

Techniques  for  resolving  the  sign  of  the  line-of-sight  velocity  component  are 
available.  They  are  not  discussed  here  because  they  were  not  available  for 
the  subject  test. 
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aerosol  provides  sufficient  scattering,  thus  enabling  operation  in  "clear  air" 
conditions;  and  (4)  the  ambient  aerosol  tracer  has  a small  inertia  and  re- 
sponds quickly  to  variations  in  airspeed  and,  thus,  can  be  a good  turbulence 
indicator. 


A useful  instrument  must  also  incorporate  means  for  scanning  the 
focal  volume  through  a selected  sequence  of  points  in  space  and  to  effect 
the  required  signal  processing,  on-line  readout,  and  permanent  recording 
requirements.  The  hardware  implementation  of  the  mobile  laser  Doppler 
unit  used  during  this  investigation  is  discussed  in  the  following  sub- 
sections. The  overall  configuration  is  summarized  in  Fig.  2-1. 

2.1.1  Basic  Laser  Doppler  Velocimeter  Optical  System 

The  optical  system  is  of  monostatic  design  and  utilizes  a continuous- 
wave  laser.  The  arrangement  depends  on  focusing  the  transmitter  telescope 
at  the  location  of  interest  for  its  spatial  resolution  property.  Details  of  the 
optical  arrangement  are  shown  in  Fig.  2-2. 


A horizontally  polarized,  20-watt,  continuous -wave  CO^  laser  beam 
(10.6-micron  wavelength)  emerges  from  the  laser  (T)  and  is  deflected  90 
degrees  by  a mirror  (T)  and  by  a 90%  reflecting  beamsplitter  (T).  The  ap- 
proximately 6-mm  diameter  beam  then  passes  through  a Brewster  window  (7) 
and  a CdS  quarter  waveplate  (5)  which  converts  the  beam  to  circular  polari- 
zation. The  beam  impinges  on  the  secondary  mirror  (fT)  and  is  expanded  and 
reflected  into  the  primary  mirror  (30-cm  diameter)  (T)and  then  focused  out 
into  the  atmosphere.  A small  portion  of  the  original  laser  beam  is  transmitted 
through  the  beamsplitter  (T)  and  is  used  as  a local  oscillator  after  being  ro- 
tated to  vertical  polarization  by  a half  waveplate  (9).  Energy  scattered  by 
aerosols  at  the  focal  volume  (T5)  is  collected  by  the  primary  mirror  Q),  col- 
limated by  the  secondary  (^6),  and  passed  through  the  quarter  waveplate  (li). 

A wire  stop  (l6)  eliminates  most  of  the  secondary  mirror  reflection  of  the  out- 
going beam.  The  quarter  waveplate  changes  the  polarization  of  the  aerosol  back- 
scattered  radiation  from  circular  to  vertical  linear  polarization.  The  vertically 
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FIGURE  2-1.  SYSTEM  CONFIGURATION. 


FIGURE  2-2.  TYPICAL  OPTICAL  COMPONENT  CONFIGURATION  OF  LOCKHEED 
LASER  DOPPLER  VELOCIMETER. 
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polarized  beam  is  approximately  78%  reflected  off  the  Brewster  window  (?)  and 
directed  through  the  beamsplitter  (lo)  where  it  is  combined  with  the  local  oscil- 
lator radiation.  After  passing  through  the  collecting  lens  (0),  the  two  beams 
are  photomixed  on  the  detector  ^2)  in  a heterodyne  configuration.  The  electrical 
output  of  the  detector  (12)  is  amplified  (14)  with  a 5-MHz  bandwidth,  20-dB  gain 
low-noise-type  preamplifier  and  fed  into  a spectrum  analyzer  (l5) . 


An  alternative  operating  configuration  consisted  of  utilizing  the  portion 
of  the  outgoing  beam  backscatter ed  into  the  interferometer  by  the  secondary 
mirror  (?)  as  the  local  oscillator  beam.  This  mode  of  operation  is  less  sus- 
ceptible to  optical  misalignment  difficulties  and  was  the  technique  used 
during  this  investigation.  When  incorporated,  the  optical  leg  (F)(2)  (9}  was 
deactivated  by  the  beam  stop  (^  and  the  wire  stop  (l^  was  removed. 

2.1.2  Optical  Scanning  System 

In  order  to  provide  the  flexibility  to  operate  in  various  modes  for  which 
the  system  was  designed,  a scanning  arrangement  as  shown  in  Fig.  2-3  is  used. 
Modes  of  operation  include  vortex  tracking  (not  required  for  the  measurements 
described  in  this  document)  and  Velocity  Azimuth  Display  (VAD)  for  the  meas- 
urement of  atmospheric  wind.  The  mirror  assembly,  AB,  can  be  rotated  about 
the  vertical  axis  for  the  scanning  in  azimuth  necessary  for  the  VAD  (also  called 
the  conical  scan  mode  of  operation).  Mirror  A is  adjusted  to  control  the 
elevation  angle  of  the  beam,  thus  controlling  the  cone  angle  of  the  conical 
scan.  The  scanning  hardware  as  deployed  on  the  mobile  van  is  shown  in 
Fig. 2-4. 

Range  scanning  of  the  system's  focal  volume  is  accomplished  by  varying 
the  distance  between  the  telescope  secondary  mirror,  E,  and  the  primary  mir- 
ror, D.  This  is  effected  by  varying  the  position  of  the  mirror,  E,  in  a con- 
trolled manner  by  an  electric-motor/optical-encoder  combination. 

The  operator  inputs  for  the  scanning  system  are  made  through  a con- 
trol panel  incorporating  thumbswitch  controls  and  light-emitting-diode 
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monitors.  The  system's  scan  capabilities  are  summarized  in  Fig.  2-5.  A 
view  of  the  mobile  unit  is  presented  in  Fig.  2-6.  Interior  views  of  the  mobile 
unit  are  shown  in  Figs.  2-7  and  2-8. 

2.1.3  Signal  Processing  System 

The  Doppler  frequency  shift  of  the  photodetector  output  is  determined 
through  the  use  of  a sampled  spectrum  analyzer  which  provides  frequency 
spectra  (intensity  of  returned  signal  as  a function  of  Dopplci  shift)  at  a rate 
of  70  signatures  per  second.  A typical  Doppler  wind  signature  is  shown  in 
Fig.  2-9.  To  yield  a line -of-sight  velocity  estimate,  a voltage  which  has  the 
same  time  behavior  as  the  Doppler  shift,  f^,  given  by  the  peak  of  the  spectrum 
is  generated. 

The  implementation  of  this  technique  is,  in  essence,  a recursive  com- 
parison method.  The  spectrum  analyzer  scan  is  driven  by  a sawtooth  voltage 
derived  from  a D/A  converter,  the  input  to  which  is  counterclocked  at  a con- 
stant rate,  hence  the  digital  output  of  the  counter  represents  frequency  on  a 
linear  scale.  At  each  new  count,  the  spectrum  analyzer  output  is  converted 
to  a digital  representation  by  an  A/D  converter,  and  the  binary  number  repre- 
senting the  current  sample  is  compared  with  that  obtained  on  the  previous 
count.  If  the  current  sample  is  the  larger  of  the  two,  it  is  saved  by  storing 
in  a latch,  along  with  the  binary  number  representing  its  frequency:  if  it  is 
smaller,  the  previous  one  is  retained  until  the  next  comparison.  This  process 
is  continued  for  the  entire  sweep.  It  is  evident  that  the  number  remaining  in 
the  frequency-store  latch,  when  the  sweep  is  completed,  corresponds  to  the 
highest  signal  power  observed;  i.e.,  the  peak  of  the  spectrum.  At  the  end  of 
each  sweep,  the  new  peak  frequency  replaces  that  obtained  on  the  previous 
sweep.  An  example  of  the  output  is  shown  in  Fig.  2-10  for  the  case  of  an  FM 
signal  of  center  frequency  2.0  MHz  (fQ)  modulated  sinusoidally  to  +200  kHz 
about  fQ  at  a 5-Hz  rate.  The  raw  spectral  information  (output  of  the  spectrum 
analyzer)  is  also  made  available  to  the  Systems  Engineering  Laboratories  (SEL) 
8 1 OA  data-logging  minicomputer  which  is  programmed  to  generate  its  own  esti- 
mate of  the  spectral  peak. 
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FIGURE 


FIGURE  2-6.  LOCKHEED-HUNTSVILLE  MOBILE  ATMOSPHERIC  UNIT. 
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FIGURE  2-7.  INTERIOR  VIEW  OF  MOBILE  ATMOSPHERIC  UNIT  LOOKING  FORWARD. 


FIGURE  2-8.  INTERIOR  VIEW  OF  MOBILE  ATMOSPHERIC  UNIT. 
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FIGURE  2-10.  OUTPUT  OF  SIGNAL  PROCESSOR  FOR  FREQUENCY 
MODULATED  INPUT. 


A provision  is  included  for  tracking  single-sideband  suppres sed-carrier 
signals,  with  an  identification  of  upper  or  lower  sideband  such  that  when  used 
in  conjunction  with  an  acousto-optic  modulator  the  unit  can  discriminate  the 
sign  of  the  Doppler  shift.  The  signal  feedthrough  at  the  translated  frequency 
can  also  be  discriminated  against  digitally,  thus  eliminating  the  need  for  a 
notch  filter.  Such  a translating  device  was  not  used  in  the  data  collection  de- 
scribed herein. 

2.1.4  Data  Recording  and  Display 

The  primary  data-gather ing  function  is  performed  by  the  SEL  810A 
general-purpose  minicomputer.  Data  acquired  by  the  Mobile  Atmospheric 
Unit  is  formatted  by  the  computer  software  and  stored  on  magnetic  tape  for 
subsequent  off-line  processing.  The  SEL  7-track  tape  control  and  magnetic 
tape  units  allow  digital  recording  of  data  at  800  bpi  at  45  ips.  The  data 
logged  by  the  computer  include: 
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All  scan  volume  location  parameters 
"Mode  of  operation"  identifier 

The  instantaneous  line-of- sight  velocity  information 
The  Doppler  spectrum  peak  strength 
Full  spectrum  intensity  and  frequency  information 
A data-quality  identifier. 

Properties  of  the  Doppler  spectrum;  namely,  the  amplitude  and  frequency 
corresponding  to  the  spectral  peak,  are  obtained  as  a result  of  on-line  com- 
puter processing  except  for  the  frequency,  which  is  also  obtained  by  the 
spectral  peak  locator  (velocity  processor)  discussed  previously.  The  latter 
allows  some  flexibility  for  on-line  operator  displays  (see  below).  The  veloc- 
ity processor  estimate  of  the  instantaneous  line-of-sight  velocity,  updated  at 
a 70-Hz  rate,  is  available  in  analog  format  which  can  be  recorded  directly  on 
a stripchart  recorder,  an  option  which  is  extremely  useful  during  the  VAD 
mode  of  operation  for  monitoring  the  characte ristic  profile. 

A view  of  the  computer  and  associated  LDV  electronics  is  shown  in 
Fig. 2- 11. 

2.2  WINDS  ALOFT  SENSING 

Using  the  basic  system  outlined  previously  it  is  possible,  by  scanning 
operations,  to  determine  the  three-component  wind  field  at  any  altitude  be- 
tween 16  and  865  meters.  The  scanning  method  employed  is  commonly  re- 
ferred to  as  the  Velocity  Azimuth  Display  (VAD)  technique  and  was  first  used 
by  I.hermitte  and  Atlas  in  conjunction  with  a microwave  radar  (Ref  , 6). 

The  telescope  is  focused  at  the  altitude  of  interest,  the  beam  being  di- 
rected at  a zenith  angle,  a.  The  beam  is  then  scanned  in  azimuth,  thus  tracing 
out  a circle  at  the  selected  altitude  (Fig.  2-12). 


Amplitude  Threshold 
Setting 


Frequency  Threshold 
Setting 


SEL  810A 
Computer 


Teletype 


Spectrum 
Analy.  er 


O s c il  1 o s c 


FIGURE  2-11.  COMPUTER  MAINFRAME  TELETYPE  AND  LASER 
DOPPLER  VELOCIMETER  ELECTRONICS. 
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FIGURE  2-12.  PRINCIPLE  OF  VELOCITY  AZIMUTH  DISPLAY  OPERATION. 
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The  instantaneous  line-of-sight  component  of  velocity  within  the  sensing 
volume  as  measured  by  the  LDV,  v , is  given  by 

v = v,  sina  cos  (0-0  ) + w cosa  , 
r h o 

v^  and  ©o,  respectively,  being  the  speed  and  direction  of  the  horizontal  wind 
component,  and  w the  vertical  component  at  the  height  being  sampled.  The 
azimuthal  dependence  of  v^  is  sufficient  to  yield  the  horizontal  speed  and 
direction  and  the  vertical  component  of  velocity. 

In  the  present  mode  of  operation,  the  system  is  unable  to  distinguish 
between  positive  and  negative  values  of  vr-  Therefore  it  is  the  absolute 
value  of  vr  (|vr|)  that  is  sensed.  This  results  in  a signal  as  shown  in  Fig. 
2-13  instead  of  the  sinusoidal  signal  as  shown  in  Fig.  2-12.  Thus,  there  is 
an  ambiguity  of  180  deg  in  the  wind  direction  since  it  is  uncertain  which  peak 
in  Fig.  2-13  represents  looking  into  the  wind.  In  practice,  the  ambiguity  is 
removed  by  the  operator  recording  approximate  wind  direction.  The  data 
processing  technique  can  then  calculate  exact  wind  direction.  This  resolves 
all  wind  direction  ambiguities  if  the  operator's  input  estimate  is  within  +89 
deg  of  the  true  wind  direction. 
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FIGURE  2-13.  AZIMUTH  ANGLE  DEPENDENCE  OF  MEASURED 
VELOCITY  COMPONENT 
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While  operating  in  the  VAD  mode,  the  system  is  capable  of  measuring 
winds  at  n (n  = 1 through  8)  altitudes  (that  can  be  dialed  in  by  using  thumb- 
switches)  in  sequence  over  a total  time  period  of  5np  seconds,  where 

5 sec  = time  for  conical  sweep  at  one  altitude  for  this  test, 

n = number  of  altitudes  to  be  interrogated,  and 

p = number  of  VAD  scans  at  each  altitude 

(can  be  chosen  to  be  1 through  7). 

During  this  investigation,  n = 8,  p = 1 were  utilized,  thus  allowing  the  meas- 
urement of  wind  at  eight  altitudes  every  40  seconds.  Sequencing  through  the 
eight  altitudes  was  performed  for  the  test  duration  required.  Data  were  also 
taken  with  continuous  scanning  at  one  altitude  (n  = l,  p=l). 
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COMPUTER-SOFTWARE  SYSTEM  DEVELOPMENT 


Acquisition  and  processing  of  the  LDV  signature  is  accomplished  by 
means  of  a compact  data-handling  system  developed  specifically  for  the  Lock- 
heed-Hunts ville  MAU.  The  general  elements  of  the  MAU  data  acquisiton-and- 
processing  system  are  shown  in  Fig.  3-1.  The  digitized  LDV  intensity 
versus  frequency  signal  along  with  its  coordinates  in  space  is  fed  into  the 
SEL  810  minicomputer.  Preprocessing  of  the  LDV  signal  is  carried  out  on 
the  minicomputer  utilizing  on-line  computer  programs  written  in  SEL  ma- 
chine language.  Information  from  the  SEL  810  is  stored  on  magnetic  tape 
and  is  used  as  an  input  to  the  off-line  processing  algorithms.  Off-line  proc- 
essing of  the  LDV  signal  is  carried  out  on  a Univac  1108  computer  with 
programs  written  in  FORTRAN  language  and  using  card  inputs  with  informa- 
tion from  the  data  logs  to  supplement  the  data.  The  final  output  consists  of 
printouts  and  plots.  A description  of  the  data  logger  program  and  the  VAD 
program  and  the  operational  characteristics  of  these  programs  is  given  in 
the  following  sections. 

3.1  DESCRIPTION  OF  LASER  DOPPLER  VELOCIMETER  SOFTWARE  SYSTEM 

Data  acquisition  in  the  MAU  is  carried  out  by  the  SEL  Data  Logger 

program.  The  Data  Logger  program  preprocesses  and  records  the  LDV 

signal.  A flowchart  of  the  Data  Logger  program  is  given  in  Fig.  3-2.  A 

sweeping  spectrum  analyzer  is  used  to  detect  the  Doppler  shift  frequency. 

A diagram  of  the  output  of  the  spectrum  analyzer  is  shown  in  Fig.  3-3.  For 

each  sweep  of  10-,  20-,  or  50-millisecond  duration,  the  Data  Logger  saves  the 

maximum  amplitude  LDV  signal,  I , and  its  corresponding  frequency,  Vmg , 

which  are  above  both  the  amplitude  and  frequency  thresholds.  The  definition 

of  I and  V and  the  shape  of  the  characteristic  LDV  spectrum  are  shown 
ms  ms  r 

in  Fig.  3-3.  The  velocity  at  maximum  signal  intensity  is  vms  and  is  taken 


FIGURE  3-2.  DATA  LOGGER  MACRO  FLOWCHART. 
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FIGURE  3-3.  TYPICAL  SPECTRUM  ANALYZER  OUTPUT  IN  VELOCITY 
AZIMUTH  DISPLAY  SCAN. 


to  be  the  characteristic  line-of-sight  velocity  associated  with  the  flow  phe- 
nomenon. This  choice  is  made  because  the  maximum  signal  intensity  is 
obtained  from  the  center  of  the  focal  volume. 

i 

The  output  from  the  Data  Logger  program  consists  of  as  a func- 

tion of  time  and  space.  From  the  output  of  the  Data  Logger,  the  wind  field 
can  be  reconstructed  using  off -line  processing  routines. 

Final  processing  of  the  LDV  measurements  is  carried  out  by  the  VAD 
program.  A macio  flowchart  of  the  VAD  program  is  shown  in  Fig.  3-4.  In 
this  off-line  program  the  array  of  V^s  values  which  is  a function  of  time 
and  space  is  processed  to  yield  the  three-dimensional  wind  field  in  the  VAD 
mode.  The  program  is  formulated  to  calculate  winds  for  both  the  translated 
and  non-translated  LDV  signal.  A translated  signal  is  provided  when  the  LDV 
system  includes  a frequency  translator  which  distinguishes  between  positive 
and  negative  values  of  line-of-sight  velocity.  A non-translated  signal  provides 
only  the  absolute  value  of  line-of-sight  velocity  as  shown  in  Fig.  2-13.  How- 
ever, during  the  course  of  this  research  effort,  all  of  the  data  acquisition- 
and-processing  were  done  in  the  non-translate  mode. 

Three  techniques  have  been  implemented  to  compute  the  three- 
dimensional  wind  components:  (1)  a peak  algorithm  where  the  magnitude  and 
location  of  the  peak  signal  in  the  sinusoidal  LDV  VAD  signature  are  used  to 
compute  the  velocity  components;  (2)  a spectral  processing  for  the  winds 
using  the  derectified  signal;  and  (3)  a sine  curve  fit.  The  final  output  is  a 
printout  (and  selected  plots)  of  the  wind  velocity  components  as  a function 
of  altitude  and  time.  Wind  velocity  components  are  given  for  both  rectilinear 
orthogonal  components  and  cylindrical  components. 

3.2  OPERATION  OF  LASER  DOPPLER  VELOCIMETER  SOFTWARE  SYSTEM 

Operation  of  the  Lockheed-Huntsville  MAU  involves  initialization  of  the 
SEL  Data  Logger  program  and  the  recording  of  the  LDV  signatures.  After 
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the  measurements  have  been  recorded  at  the  proper  threshold  settings,  the 
VAD  program  is  used  to  process  the  data  into  VAD  plots.  It  is  useful  to  ex- 
amine the  data-processing  operations  involved  in  determining  the  three- 
dimensional  wind  velocity  field  from  the  LDV  measurements  in  terms  of 
on-line  preprocessing  and  off-line  processing. 

3.2.1  On-Line  Data  Processing 

The  conical  scan  VAD  measurements  are  preprocessed  on  the  SEL 
computer  and  recorded  on  magnetic  tape.  The  output  from  the  SEL  computer 
consists  of  the  basic  V^s  signal  as  well  as  additional  test  parameters  which 
are  listed  in  Table  3-1.  It  is  useful  to  consider  the  type  of  data  recorded  by 
the  SEL  during  normal  operations  in  the  VAD  mode. 


A dump  of  a sample  output  tape  from  the  Data  Logger  operating  in  the 
VAD  mode  is  shown  in  Fig.  3-5,  each  row  of  data  corresponding  to  information 
recorded  for  each  spectrum-analyzer  sweep.  The  information,  packed  into 
four  computer  words  on  tape,  is  separated  into  13  columns  in  the  printout  with 
interpretation  of  the  various  columns  being  as  follows: 


Column 


Interpretation 


1 (PEAK/MAX)  An  indication  as  to  whether  the  peak  Doppler  frequency 

(fpk  corresponding  to  Vpk  in  Fig.  3-4)  or  maximum 
signal  Doppler  frequency  (fms  corresponding  to  Vms  in 
Fig.  3-4)  is  used.  Indication  is  0 if  maximum  signal 
frequency  is  used  and  is  1 if  peak  velocity  signal  is 
used. 

2 (0  = 0)  A conical-scan  azimuth  reference  which  is  nonzero  when 

the  reference  switch  is  activated. 

3 (NOR/VAD)  An  indication  as  to  whether  the  system  is  operating  in 

the  normal  (vortex  track)  (1)  or  VAD  (0)  mode. 

4 (SWEEP  SPEED)  Sweep  speed  of  the  spectrum  analyzer  trace  in  msec/cm 

5 (DIG  TRCK)  The  computer-calculated  estimate  of  f , (percent  of 

full  scale).  p 

6 (DATA  ACCEPT)  An  indication  that  there  is  (1)  or  is  not  (0)  an  output 

above  the  frequency  and  amplitude  thresholds  during  a 
sweep. 
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TABLE  3-1 

INPUT  FROM  SYSTEMS  ENGINEERING  LABORATORIES  COMPUTER 


Each  VAD  run  is  in  a separate  file.  The  following  data  are  recorded 
for  each  spectrum  analyzer  sweep. 

1. 

Spectral  Sample  Count  Across  Spectrum  Analyzer  S /eep, 

Corresponding  to  V 

r ° max 

2. 

Amplitude  at  Above  Point 

3. 

Data  Acceptable  Flag 

4. 

Flag  for  Spectrum  Analyzer  Sweep  Speed 

5. 

Flag  for  Translator 

6. 

Flag  for  Positive  or  Negative  Frequency  (used  only  when 
translator  is  used) 

7. 

Flag  for  V .or  V 
b pk  ms 

8. 

Flag  for  Conical  Scan  or  Normal  Scan 

For  Conical  Scan 

9. 

Height  above  Van 

10. 

Flag  for  Azimuth  Switch 

11. 

Cone  Angle 

k. 


Column 


Interpretation 


/ 


7 ( + /-) 


8 (TRANS/NON- 
TRANS) 


9 (SPECTRUM 
INTENSITY) 


An  indication  as  to  the  sense  of  the  Doppler  shiiL: 

1 = moving  toward  MAU,  0 = moving  away  from  MAU. 


An  indication  as  to  whether  or  not  a frequency  transla- 
tor was  incorporated.  During  this  investigation,  it  w.ii 
not  incorporated.  1 = Yes,  0 = No. 

Peak  amplitude  of  the  Doppler  spectrum  in  region  above 
a frequency  threshold  (arbitrary  units  ranging  from  0 
to  1024). 


10  (NBR  ROTN)  Number  of  successive  VAD  scans  for  a particular  altitude. 

11  (ALT  (M))  Altitude  of  VAD  for  particular  sweep  in  meters. 

12  (LTRNC  TRCK)  On-line  frequency-tracker  estimate  of  f , (should  be 

approximately  equal  to  column  5). 


13  (CONE  ANGLE)  Half-angle  of  VAD  cone  in  degrees. 


3.2.2  Off-Line  Data  Processing 

Three  algorithms  for  calculating  the  mean  wind  speed  and  direction  from 
the  VAD  signature  have  been  developed.  For  each  of  these  algorithms,  mean 
wind  speed  and  direction  are  calculated  for  each  5-second  VAD  sweep.  Stan- 
dard deviations  of  wind  speed  and  direction  can  be  calculated  from  multiple 
VAD  sweeps.  The  data  output  of  the  LDV  system  operating  in  the  VAD  mode 
are  line-of-sight  velocities  measured  at  a selected  number  (usually  350  in 
the  current  Lockheed  system)  of  distinct  points  around  the  VAD  cone.  Recall 
that  the  line-of-sight  velocity  signature  is  theoretically  sinusoidal  in  the  VAD 
mode  (cf.,  Fig.  2-12). 

For  all  the  processing  algotithms,  preprocessing  of  the  data  occurs; 
the  preprocessing  includes: 

a.  Save  line-of-sight  velocities  for  one  roation  of  scanner. 

b.  If  two  or  more  rotations  occur  at  the  same  altitude,  average 
with  previous  rotation(s). 

c . Assign  azimuth  angle  to  each  point  (assuming  constant  rotation 
rate). 
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d.  Edit  data  to  eliminate  spurious  points.  Eliminated  points  are  de- 
clared unacceptable. 

e.  Plot  line-of-sight  velocity  versus  azimuth  angle. 


The  edit  criterion  for  the  elimination  of  spurious  points  is  that  the  i 
point  is  eliminated  if 


v -v  ...  >.2v  .... 

r,  1 r, l + l r , l+l  ’ 


v .-v  . , > .2  v ... 

r,i  r,i-l  r,  i-  1 


A sample  plot  of  unedited  line-of-sight  velocity  versus  azimuth  angle  is  shown 
in  Fig.  3-6.  Points  for  unacceptable  data  are  shown  as  zero  velocity. 


3.2.2. 1 Peak  Algorithm 


For  the  calculation  of  wind  velocity  by  the  peak  algorithm,  the  pro- 


cedure is: 


a.  Filter  data  with  an  n-point  moving  average. 


b.  Identify  the  two  peak  velocity  points,  V j and  V that 
at  a minimum  of  90  degrees  apart.  P P 


c.  Compute  the  magnitude  of  the  horizontal  component  of  the  wind 


Vh  = 


V . + V , 
pi  p2 


2 sina 


d.  Compute  horizontal  wind  angle  with  help  of  estimated  wind 
direction. 


e.  Compute  the  vertical  component  of  the  wind 


w = 
v 


V , - V , 
pi  p2 


2 cosa 


f.  Derectify  VAD  signal  if  no  translator  is  present  and  plot  derectified 
signal.  The  signal  is  derectified,  so  that  the  positive  peak  of  the  derectified 
signature  is  the  peak  which  is  closer  to  the  estimated  wind  direction  recorded 
by  the  operator  at  the  time  the  data  are  measured. 


The  purpose  of  the  filter  is  the  elimination  of  nonuniformities  in  the 
data.  This  is  especially  important  for  the  peak  algorithm.  The  individual 
data  points  are  filtered  with  an  n-point  moving  average  filter  (n  is  usually 
21).  Thus, each  line-of-sight  point  is  filtered  by 


V 

r,  i 


V 

. r. 


n-  1 
2 


+ 


. + V 

r, 


where  ^ is  the  filtered  value  of  the  line-of-sight  velocity  to  be  used  in 
further  calculations.  A plot  of  the  filtered  line-of-sight  velocities  for  a 21- 
point  filter  is  shown  in  Fig.  3-7.  Additional  samples  of  LDV  signatures  (includ- 
ing raw  data,  filtered  data,  and  derectified  data)  are  presented  in  Appendix  A. 

When  the  LDV  data  are  measured,  an  approximation  of  the  wind  angle 
is  recorded.  The  calculated  wind  angle  is  the  azimuth  angle  of  the  peak 
which  is  closer  to  the  estimated  wind  angle.  For  small  values  of  vertical 
component  of  velocity  the  wind  angle  plus  90  deg  is  the  angle  at  which  the 
line-of-sight  velocity  is  theoretically  zero.  This  angle  is  used  for  the  de- 
rectification of  the  line-of-sight  signal.  A plot  of  the  derectified  (edited  but 
unfiltered)  line-of-sight  velocity  is  presented  in  Fig.  3-8. 


3. 2. 2. 2 Fourier  Coefficient  Algorithm 

The  Fourier  coefficient  (or  spectral)  algorithm  computes  the  funda- 
mental harmonic  of  the  line-of-sight  velocity.  The  Fourier  series  for  a 
generalized  periodic  function  is 


cx  - , 

f(x)  = -=-  + > (a  cosnx  + b sinnx), 

c.  n n 

n = 1 

where  the  Fourier  coefficients  are  given  by 

L 

2 V"*  \r  ,2?rin. 
an  : L L Vi  C°s(-L_ )■ 
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FIGURE  3-7.  FILTERED  LINE-OF-SIGHT  VELOCITY  FOR  VELOCITY 
AZIMUTH  DISPLAY  MODE. 
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where  the  L discrete  points  are  uniformly  spaced  between  0 and  2?r,  and  V. 

th  ^ 

is  the  line-of-sight  velocity  for  the  i point  in  the  scan.  For  the  calculation 

of  the  three  wind  components,  it  is  necessary  to  calculate  only  a^,  aj,  and  b^. 
However,  the  second  and  third  harmonics  (i.e.,  n = 2 and  n = 3)  are  also  cal- 
culated as  an  indication  of  the  nonuniformity  of  the  wind.  A correction  factor 
must  be  applied  to  compensate  for  the  absence  of  data  points  which  lie  below 
the  yelocity  threshold.  Although  the  line-of-sight  signature  in  the  VAD  mode 
is  theoretically  sinusoidal,  the  existence  of  the  velocity  threshold  causes  the 
signature  to  appear  as  shown  in  Fig.  3-9.  For  the  purposes  of  deriving  the 
correction  factor,  the  origin  is  chosen  so  that 


v 


r 


12 

2 


+ b sin9. 


In  Fourier  analysis,  a^/2  is  the  d-c  component.  If  the  vertical  component  of 
velocity  is  small  compared  with  the  magnitudeof  the  velocity  threshold,  the 
range  of  azimuth  angles  for  which  the  signature  is  zero  is  as  shown  in  Fig. 
3-9.  The  parameter,  Z,  is  the  ratio  of  the  number  of  zero  points  to  the  total 


FIGURE  3-9.  LINE-OF-SIGHT  VELOCITY  SIGNATURE  OF  VELOCITY 
AZIMUTH  DISPLAY. 
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number  of  points  in  the  VAD  scan.  Let  a*  and  b^'  be  the  values  of  the  Fou- 
rier coefficients  calculated  from  the  line-of-sight  velocity  signature.  The 
integral  form  of  the  equations  for  the  Fourier  coefficients  is  used.  For  the 
theoretical  velocity  signature, 

2.n 

ao  = lA 

0 

However,  for  the  actual  velocity  signature, 

■tr-Zit/Z  2ir-Zjr/2 

J 

^t/2 


^7 r 

a*Q  = l/j r j (aQ/2  + bj  sin9)  d9  = aQ. 


a0  = 1/t 


(a0/2  + bj  sin9)  d9  + 


/ 

t+Zt/2 


(a^/2  + bj  sin9)  d9 


= aod-Z). 


Therefore,  a correction  factor  of  l/(l-Z)  must  be  applied  to  the  calculated  value 
of  aQ. 


For  the  sine  term,  the  theoretical  velocity  signature, 

Zir 

b j = 1/7 r J"  (aQ/2  + bj  sin9)  sin9  d9  = bj- 

0 

For  the  actual  velocity  signature, 

Zn-Zn/Z 


S = iA 


■k-Ztt/Z 


J"  (a^/2  + b j sin9)  sin9  d9  + J"  (a^/Z+b^  sin9)  sin9  d9 


[Zn/Z 


7T+Z  Tt/Z 


Therefore,  a correction  factor  l/[l  - Z + (sinZir)/ir]  must  be  applied  to  the 
horizontal  wind  components. 

The  calculation  procedure  for  the  spectral  algorithm  is: 

a.  Compute  the  Fourier  coefficients  for  the  d-c  level  and  the  first 
three  harmonics  (n=l,2,3): 
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an  = 2/L  ^ V.  cos(2irin/L), 
i=  1 

L 

b = 2/L  ^ V.  sin(2jr  i n/L) . 
i=l 


b.  Compute  the  corrected  fundamental  harmonic: 

Corrected  Value  = — ta-  ■ \ — — — r-7—  * Calculated  Fundamental. 
1 - Z + (sinZjr  )/ir 

c.  Compute  the  vertical  wind  correction: 

Corrected  Value  = ^ 1 z * Calculated  Value. 

d.  Compute  the  magnitude  of  the  horizontal  component  of  velocity: 


Vh  = 


sma 


e.  Compute  the  horizontal  angle: 


Angleh  = Atanfb  ^/a  j). 


f.  Compute  the  vertical  component  of  wind  velocity: 

-ao 


2 cosa. 


3. 2. 2. 3 Sine -Cur ve -F it  Technique 

The  line-of- sight  velocity  signature  is  sinusoidal  for  a uniform  wind. 
Therefore,  the  best  sinusoidal  wave  which  fits  the  data  in  a least- 
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squares  sense  is  determined  by  finding  the  values  of  the  coefficients  A,B, 
and  C which  minimize  (summing  over  all  acceptable  data  points) 

^ ' (V.  - C - A cosfh  - B sin8.)^, 
i 

where  is  the  line-of-sight  velocity  (derectified)  at  point  i,  and  (h  is  the 
azimuth  at  point  i.  A,B,  and  C are  obtained  from 

[]CC°s20i]A  + []CCOs9i  sin9i]B  + £Ccos0i]C  = S Vi  cosei- 

i i i 

J^^cosQ.  sinQ.jA  + j^^sin20.jB  + F^sin^C  = ^ V.  sin9i> 

i i i 

[j^cose^A+^sine.jB  + nC  = £v.. 
i i i 

The  steps  for  calculating  the  wind  using  the  least-squares  algorithm  are: 

a.  Find  the  least-squares  curve  fit  for  a sine  wave  to  minimize 

Y]  (V.  - C - A cos0  . - B sin6.  )2, 
i 

where 

V.  is  line-of-sight  velocity  at  point  i, 

0|  is  azimuth  at  point  i,  and 

C,  A,  and  B are  coefficients  to  be  calculated. 

b.  Compute  the  magnitude  of  the  horizontal  component  of  velocity: 


Vh  = 


/a2  + B2 


c.  Compute  the  horizontal  angle: 


Angle^  = Atan(B/A). 
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4.  DATA  COLLECTION  AT  TABLE  MOUNTAIN 
4.1  TEST  DESCRIPTION 

During  February  and  March  1976,  the  mobile  LDV  was  deployed  at  the 
National  Oceanic  and  Atmospheric  Administration's  Table  Mountain  Test  Site, 
about  16  km  north  of  Boulder,  Colorado.  During  this  period,  the  wind  varied 
from  0 to  15  m/sec  at  the  30-m  altitude  Above  Ground  Level  (AGL)  to  0 to  18 
m/sec  at  148-m  altitude.  Light  to  moderate  snowfall  occurred  twice,  and  no 
rainfall  occurred  during  this  period  of  time.  Data  were  collected  under  atmos- 
pheric conditions  from  nearly  laminar,  thermally  stable  flow  to  strong,  gusty, 
thermally  unstable  flow.  A diagram  of  the  test  site  layout  is  shown  in  Fig.  4-1, 
and  a photograph  of  the  test  site  is  shown  in  Fig.  4-2. 

By  standard  meteorological  convention,  the  wind  direction  is  the  di- 
rection from  which  the  wind  comes.  Therefore,  when  wind  is  given  in  speed- 
direction  coordinates,  the  direction  is  the  direction  from  which  the  wind  comes, 
measured  from  north.  The  coordinates  chosen  for  the  test  are  shown  in  Fig. 
4-1.  For  wind  expressed  in  rectangular  coordinates,  the  positive  u component 
represents  wind  coming  from  the  positive  x direction  (i.e.,  from  east  to  west). 
Similarly,  a positive  v component  represents  wind  coming  from  the  positive 
y component  (i.e.,  from  north  to  south).  For  example,  a wind  coming  from 
1 35-deg  has  a positive  u component  and  a negative  v component. 

The  tower  was  instrumented  with  propeller-vane  anemometers,  Model 
8002  manufactured  by  R.M.  Young  Co.  at  approximately  30,  60,  90,  120,  and 
150  m.  The  true  altitudes  AGL  of  the  tower  instrumentation  were  31,  67,  95,  125, 
and  148  m and,  in  this  report,  are  used  interchangeably  with  the  approximate 
values.  As  reported  in  Ref.  7,  the  anemometers  were  mounted  on  booms  ex- 
tended 2.5  m north  of  the  tower  structure  and  provided  sensor  exposu:  e with 
undisturbed  airflow  except  for  a 90-deg  sector  centered  roughly  on  1 75-deg 
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FIGURE  4-1.  NATIONAL  OCEANIC  AND  ATMOSPHERIC  ADMINISTRATION 
TABLE  MOUNTAIN  TEST  SITE  LAYOUT. 
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azimuth.  A careful  calibration  was  employed  to  certify  the  tower  instrumen- 
tation measurements  prior  to  their  mounting  on  the  tower.  The  calibration 
showed  no  observable  systematic  bias  errors  and  less  than  _+10  cm/sec 
difference  between  anemometers. 

The  anemometer  outputs  were  transmitted  to  a computer-controlled 
data-acquisition  system  where  each  output  was  sampled  once  per  second  and 
the  data  processed  in  real  time  to  yield  the  time-averaged  wind  data. 

A summary  of  runs  for  which  measurements  were  made  is  shown  in 
Table  4-  1 . 

4.2  TEST  OBJECTIVES  AND  OPERATING  PROCEDURE 

The  primary  objectives  of  the  test  were  to  establish  the  validity  of  the 
VAD  technique  for  three-dimensional  wind  measurement  and  to  optimize 
data-processing  algorithms  to  improve  the  accuracy  of  VAD  wind  measure- 
ment. A secondary  objective  was  to  investigate  the  effect  of  averaging  time 
on  the  accuracy  of  VAD  wind  measurement  relative  to  the  tower  measure- 
ments. 

As  discussed  in  Section  2,  the  MAU  VAD  scan  capability  allows  wind 
measurement  to  be  made  at  eight  altitudes,  each  requiring  5 seconds  for  a 
360-degree  azimuth  rotation.  This  requires  40  seconds  for  a complete  scan- 
ning of  the  eight  altitudes.  The  test  was  initiated  with  the  laser  scanning  at 
30,  45,  60,  75,  90,  120,  150,  and  300-m  altitudes. 

4.3  DATA  ANALYSIS  OF  INITIAL  RUNS 

The  winds  calculated  with  the  three  algorithms  described  in  Section  3 
are  shown  in  Figs.  4-3  through  4-7  for  one  16-minute  sample.  The  sole 
difference  between  the  three  parts  of  each  figure  is  the  data-processing  algo- 
rithm used.  Some  general  observations  are  that  the  peak  algorithm  tends  to 
give  mean  wind  values  that  are  biased  slightly  high  compared  with  the  tower- 
measured  mean  winds.  This  is  primarily  due  to  temporal  variations  in  the 
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wind  at  higher  frequencies  than  the  fundamental  frequency  in  the  VAD  signa- 


ture; thus,  the  peak  of  the  signature  (i.e.,  the  fundamental  frequency  with 
superimposed  higher-frequency  variations)  has  a value  that  is  higher  than 
the  averaged  values.  The  interpretation  is  confirmed  by  the  fact  that  the 
runs  for  which  the  greatest  deviation  above  the  diagonal  line  occurred  are 
the  runs  for  which  the  greatest  value  of  wind  standard  deviation  (measured 
by  the  tower)  occurred.  In  particular,  for  the  points  of  Run  9 farthest  off 
the  45-degree  line  in  Figs.  4-5  and  4-6,  the  standard  deviation  of  the  wind 
speed  was  1.82  m/sec  and  1.75  m/sec,  respectively,  whereas  the  standard 
deviation  of  the  wind  speed  was  usually  between  0.9  m/sec  and  1.4  m/sec. 

It  is  also  observed  that  the  spectral  algorithm  tends  to  give  mean  wind 
values  that  are  biased  slightly  low  compared  with  the  tower-measured  mean 
winds.  It  was  observed  during  the  data  processing  that  the  magnitude  of  the 
bias  is  a function  of  the  number  of  valid  line-of-sight  velocities  obtained 
during  the  VAD  scan.  The  bias  was  insignificant  for  data  when  more  than 
50  valid  data  points  were  obtained  during  each  VAD  scan.  It  is  noted  that  a 
number  of  data  points  below  this  value  was  obtained  only  rarely  at  Table 
Mountain  and  that  the  low  number  of  data  points  should  never  occur  in  air 
with  a larger  quantity  of  natural  aerosol  than  the  very  clean  air  above  Table 
Mountain.  The  least- squares  sine  algorithm  seems  to  present  no  systematic 
bias,  and  the  scatter  of  the  data  about  the  45-degree  line  is  less  for  the  sine 
algorithm  than  for  the  other  two  processing  techniques.  A comparison  of  the 
peak  algorithm  and  the  spectral  algorithm  with  the  sine  algorithm  is  shown  in 
Figs.  4-8  and  4-9.  For  the  peak  algorithm,  the  deviation  from  the  45-degree 
line  increases  as  turbulence  increases,  and,  for  the  spectral  algorithm,  the 
deviation  from  the  45-degree  line  increases  as  the  number  of  valid  points  in 
VAD  scan  decreases  below  50. 

4.4  INITIAL  DATA  COMPARISON 

Laser  and  tower  data  shown  in  Figs.  4-3  through  4-7  consist  of  a direct 
comparison  of  16-minute  mean  wind  speed.  The  proximity  of  individual  data 
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FIGURE  4-9.  COMPARISON  OF  DATA  PROCESSING  ALGORITHMS 
AT  150-M  ALTITUDE  (16-min  Means  of  Runs  6,7,  and  9). 


points  to  the  45-degree  line  is  a direct  measure  of  the  difference  between 
the  independent  measurements. 


An  important  observation  is  that  the  scatter  in  the  data  is  somewhat 
greater  than  expected  for  all  altitudes  and  all  data-processing  algorithms. 

This  was  determined  to  be  a statistical  problem  rather  than  inherent  prob- 
lems associated  with  the  VAD  assumptions.  The  rotation  rate  of  the  azimuth 
scan  has  a period  of  5 seconds,and  the  discrete  range  scan  is  hardwired  for 
performing  scans  at  a maximum  of  eight  preselected  altitudes,  resulting  in 
40  seconds  for  a profile  to  be  mapped  out  with  one  wind  calculation  performed 
per  altitude  during  each  40-second  period.  In  a 16-minute  time  frame,  24 
points  per  altitude  are  used  in  the  average  wind  calculations.  In  contrast, 
the  propeller  vanes  are  digitized  at  one  point  per  second,  which  yields  960 
points  per  16-minute  time  frame.  It  was  believed  that  the  relatively  small 
number  of  values  contributing  to  the  mean  wind  speed  when  all  eight  altitudes 
were  utilized  is  the  major  factor  in  the  scatter  of  the  laser  data.  Accordingly, 
in  later  runs,  fixed  altitude  VAD  measurements  were  implemented  to  allow 
inclusion  of  a greater  number  of  individual  VAD  wind  measurements  in  the 
16-minute  average.  It  was  believed  that  this  would  reduce  the  scatter  of  the 
comparative  data  and,  thus,  would  allow  a more  meaningful  evaluation  of  the 
VAD  concept. 


4.5  FIXED-ALTITUDE  VAD  DATA  COMPARISON 

Ideally,  the  solution  to  the  scatter  caused  by  a low  number  of  data  points 
in  a 16-minute  average  would  be  an  increased  rotation  rate.  For  the  test,  this 
solution  was  not  possible  because  the  rotation  rate  of  the  Lockheed  LDV  is 
fixed.  However,  in  order  to  quantify  the  degree  to  which  accuracy  could  be 
enhanced  by  an  increased  number  of  scans  in  a 16-minute  average,  fixed-alti- 
tude VAD  measurements  were  made.  By  making  all  scans  at  one  altitude, 
the  number  of  scans  in  each  16-minute  average  was  increased  from  24  (as  for 
the  data  discussed  in  the  previous  section)  to  192.  This  increase  was  deemed 
necessary  in  order  to  make  meaningful  conclusions  about  the  validity  of  the 
VAD  assumptions. 
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The  MAU  did  not  have  a real-time  wind-processing  capability.  To 
alleviate  the  problems  of  the  inability  to  assess  the  quality  of  data  in  the 
field,  a simplified  wind-speed  processing  algorithm  was  implemented  on  the 
computer  in  the  tower  data-pr oces sing  van  (Ref.  8).  This  allowed  post-test 
comparisons  between  the  LDV  and  tower  data  to  be  made  on  a daily  basis, 
and  proved  to  be  a valuable  asset. 

Analysis  of  the  fixed-altitude  data  showed  that  the  peak  algorithm  ex- 
hibits a slightly  high  bias  of  the  wind  speed.  This  characteristic  was  also 
noted  in  the  multiple-altitude  data.  The  bias  is  clearly  seen  in  Fig.  4- 10. 

It  is  noted  that  the  wind-speed  scale  has  been  expanded  to  display  the  differ- 
ence. As  noted  in  Section  4.3,  the  reason  for  the  bias  is  believed  to  be  tur- 
bulent fluctuations  superimposed  on  the  sinusoidal  VAD  signature  which  would 
result  from  a uniform  wind.  The  fact  that  the  bias  is  present  in  both  single- 
and  multiple -altitude  data  indicates  that  the  bias  is  not  due  to  averaging 
over  multiple  scans,  but  is  related  to  the  process  by  which  the  wind  from  each 
scan  is  calculated.  This  fact  supports  the  hypothesis  that  the  bias  is  due  to 
turbulent  fluctuations. 

The  spectral  algorithm  was  incorporated  into  the  off-line  processing 
program  to  alleviate  this  problem.  The  concept  of  the  spectral  algorithm 
was  that  the  mean  wind  would  be  indicated  by  the  first  harmonic  (in  a Fourier 
sense),  and  turbulent  fluctuations  would  be  indicated  by  higher  harmonics. 
Turbulent  fluctuations  would  not  influence  the  magnitude  of  the  fundamental 
harmonic.  The  results  are  also  shown  in  Fig. 4-10,  where  excellent  agree- 
ment between  the  LDV  and  tower  measurements  is  observed. 

The  average  wind  speed  and  direction  as  calculated  by  the  least-squares 
sine  fit  is  also  shown  in  Fig.4-10.  Maximum  deviations  between  LDV  and 
tower  measurements  are  less  than  20  cm/sec  in  wind  speed  and  2 degrees  in 
wind  direction.  The  corresponding  standard  deviation  in  wind  speed  and  di- 
rection are  depicted  in  Fig.  4- 11,  where  maximum  deviations  are  seen  to  be 
15  cm/sec  in  speed  and  7 degrees  in  direction. 
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An  interesting  record  of  data  is  shown  in  Fig.  4 -12.  During  a one  and 
one-half  hour  period  the  wind  direction  changed  approximately  70  degrees; 
the  comparisons  in  both  mean  speed  and  direction  are  excellent  with  a maxi- 
mum deviation  of  less  than  15  cm/sec  in  wind  speed  and  4 degrees  in  wind 
direction  for  the  sine  algorithm. 

Comparative  plots  are  shown  in  Figs. 4-13  through  4-15  for  laser  versus 
tower  wind  speed  and  in  Figs.4-l6  through  4-18  for  wind  direction  comparisons. 
In  general,  excellent  agreement  is  observed.  However,  there  is  still  a slightly 
high  bias  for  the  peak  algorithm,  a slightly  low  bias  for  the  spectral  algorithm 
and  no  systematic  bias  for  the  least-squares  sine  algorithm.  It  is  noted  that 
the  150-m  altitude  generally  exhibits  greater  data  scatter  than  the  lower  alti- 
tudes. Therefore,  the  scatter  at  lower  altitudes  would  be  less  than  that  shown 
in  Figs.  4- 13  through  4-18. 

4.6  EFFECT  OF  AVERAGING  TIME 

As  expected,  the  effect  of  increasing  averaging  time  is  to  decrease  the 
scatter  of  the  data  about  the  45-degree  comparison  line.  This  is  illustrated 
for  both  wind  speed  and  direction  for  the  sine  algorithm  in  Figs.  4-19  through 
4-24.  The  one-minute  averages  (i.e.,  twelve  data  points  in  the  average)  ex- 
hibit a fair  amount  of  scatter.  However,  the  six-minute  average  (i.e.,  72  data 
points  in  each  average)  shows  significantly  less  scatter.  The  15-minute 
averages  exhibit  almost  no  scatter. 

Similar  trends  were  also  observed  for  the  peak  and  spectral  algorithms. 
The  data  are  presented  in  Appendix  B in  a manner  similar  to  that  of  Figs. 4-19 
through  4-24. 

Much  of  the  data  scatter  is  attributable  to  the  distance  between  the 
meteorological  tower  and  the  Mobile  Atmospheric  Unit,  rather  than  to  in- 
accuracies in  the  measurements  made  by  the  Mobile  Atmospheric  Unit.  It 
is  expected  that  over  a long  period  of  time,  the  average  winds  at  the  two 
locations  will  be  approximately  equal.  However,  for  short-time  periods,  the 
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FIGURE  4-13.  COMPARISON  OF  LASER  AND  TOWER  15-MIN  MEAN 
WIND  SPEED  USING  PEAK  ALGORITHM. 
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FIGURE  4-17.  COMPARISON 
WIND  DIRECTION  USING  SPE 


two  measurements  would  not  be  expected  to  be  equal,  especially  when  the 
standard  deviation  of  wind  speed  is  10%  of  the  mean  wind  speed.  This  is 
illustrated  by  the  fact  that  the  data  scatter  about  the  45-degree  line  is 
greater  for  the  one -minute  means  than  for  longer  averaging  times. 

4.7  EFFECT  OF  SNOW 

The  presence  of  snow  during  LDV  wind  measurements  creates  a nega- 
tive vertical  component  of  velocity,  but  does  not  affect  the  measurement  of 
horizontal  components  of  velocity.  A typical  VAD  signature  for  light  snow 
is  shown  in  Fig.  4-25.  The  signature  is  taken  from  Run  27  for  which  the  snow- 
fall was  approximately  1 cm/hour.  There  is  a distinct  double  signature.  The 
upper  signature  is  due  to  laser  reflection  from  atmospheric  aerosol.  The 
lower  signature  is  due  to  laser  reflection  from  snow.  For  light  snow,  reflec- 
tions from  both  atmospheric  aerosol  and  snow  are  detected.  If  a snow  par- 
ticle is  in  the  focal  volume,  the  maximum  signal  intensity  is  probably  derived 
from  reflection  from  the  snow  particle.  If  there  is  no  snow  particle  in  the 
focal  volume,  the  maximum  signal  intensity  comes  from  a particle  of  the 
atmospheric  aerosol. 

All  of  the  data- processing  algorithms  average  the  two  signatures, 
weighted  by  the  number  of  points  in  each  of  the  signatures.  The  result  is  a 
signature  which  lies  between  the  "clear  air"  and  snow  signatures.  Since  the 
amplitude  and  phase  angle  of  the  sine  wave  are  unaffected  by  this  process, 
the  wind  speed  and  direction  are  unaffected.  For  example,  in  Fig.,.  4- 13  through 
4-18,  the  data  from  Run  27  exhibit  no  distinguishing  characteristics.  However, 
a negative  component  of  vertical  velocity  is  superimposed  on  the  wind  velocity. 
For  example,  the  vertical  component  of  velocity  for  Run  27  is  approximately 
-1  m/sec,  whereas  that  for  Run  23  (at  the  same  altitude  and  approximately 
equal  wind  speed),  the  vertical  component  of  velocity  is  approximately  -0.1 
m/ sec. 
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For  heavy  snow  (Runs  6 and  7),  a doable  signature  does  not  occur,  but  a 
definite  vertical  component  of  velocity  is  observed.  A typical  VAD  signature 
in  heavy  snow  is  shown  in  Fig.  4-26.  The  vertical  component  of  velocity  is 
approximately  -1  m/sec  for  Run  6. 

For  run  27,  the  variation  in  returned  laser  intensity  is  much  greater  than 
is  normally  experienced.  For  example,  a typical  example  of  a sequence  of  rela- 
tive signal  intensities  on  a linear  scale  is:  ...444,  168,  680,  272,  190,  576, 

238,  154,  190,  142,  176,  160,  304,  574,  384,  222,  160,  228,  144,  320, 

576,  352,  190,  320,  120,  144,  312,  224,  252.  . . It  is  hypothesized  that  the 

large  values  of  intensity  are  reflection  from  snow,  and  smaller  values  of  in- 
tensity are  reflection  from  natural  ambient  aerosol.  The  higher  intensities 
usually  exhibit  a line-of-sight.  velocity  which  is  discernibly  different  from  its 
adjacent  points  with  lower  intensities,  but  this  characteristic  is  not  universal. 
There  are  points  for  which  the  line-of-sight  velocity  is  discernibly  different 
from  its  adjacent  points,  but  the  signal  intensity  is  small.  It  is  possible  to 
identify  points  which  may  be  reflections  from  snow  based  upon  deviations  in 
line-of-sight  velocity  or  signal  intensity.  However,  it  is  not  possible  to  identify 
individual  points  positively  as  reflection  from  snow  on  the  basis  of  a one-to-one 
correspondence  between  deviations  in  line-of-sight  velocity  and  deviations  in 
signal  intensity. 
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5.  CONCLUSIONS  AND  RECOMMENDATIONS 

The  laser  Doppler  velocimeter  using  the  VAD  technique  is  an  accurate 
and  reliable  method  for  the  remote  measurement  of  atmospheric  wind.  Three 
data-processing  algorithms  have  been  investigated  in  order  to  determine  the 
optimal  algorithm.  It  appears  that  the  least-squares  sine-fit  algorithm  pro- 
vides the  most  consistent  data  compared  with  anemometer  data  from  a mete- 
orological tower.  The  peak  algorithm  yields  wind-speed  data  which  are 
slightly  above  those  measured  by  the  tower-mounted  anemometers  and  given 
by  the  sine  algorithm.  The  Fourier  algorithm  yields  wind-speed  data  which 
are  1 to  1.5  m/sec  below  that  given  by  the  tower  anemometers  and  the  sine 
algorithm. 

Based  on  comparisons  between  the  tower  anemometers  and  the  laser- 
measured  wind  as  calculated  by  the  least-squares  sine  algorithm  and  shown 
in  Appendix  B,  the  accuracy  of  the  LDV  varies  from  + 2 m/sec  in  wind  speed 
and  +20°  in  wind  azimuth  for  one-minute  averages  to  +1  m/sec  in  wind  speed 
and  +10°  in  wind  azimuth  for  15-minute  averages. 
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Appendix  A 

LASER  DOPPLER  VELOCIMETER  VELOCITY 
AZIMUTH  DISPLAY  SIGNATURES 

(Runs  21-28) 

This  appendix  contains  sample  velocity  azimuth 
display  signatures  to  define  data  quality  of  each 
data  run  taken  at  Table  Mountain  test  site.  The 
signatures  are  presented  in  sets  of  three  figures. 
The  first  figure  in  each  set  is  the  raw  data.  The 
second  figure  is  the  edited  and  filtered  data,  and 
the  third  figure  is  the  derectified  signature  (edited, 
but  unfiltered). 
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FIGURE  A-l.  SAMPLE  LASER  DOPPLER  VELOCIMETER  SIGNATURES. 
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FIGURE  A-l  (Continued). 
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FIGURE  A- 1 (Continued). 
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FIGURE  A-l  (Continued). 
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FIGURE  A-  1 (Concluded). 
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Appendix  B 

LASER  DOPPLER  VELOCIMETER  AND  ANEMOMETER 
WIND  COMPARISONS 

(Runs  21-28) 

The  least-squares  sine  algorithm  is  used  as  the  basis  of 
comparison  for  winds  measured  by  the  LDV  and  those  mea- 
sured by  the  tower-mounted  anemometers.  Wind-speed-and 
direction  comparisons  are  given  for  1-,  6-,  and  15-minute 
averages.  No  anemometer  data  were  available  for  Run  24. 
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FIGURE  B-l.  COMPARISON  OF  WINDS  MEASURED  BY  LASER 
DOPPLER  VELOCIMETER  AND  TOWER  ANEMOMETER. 
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FIGURE  B-l  (Continued). 
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FIGURE  B-l  (Continued). 
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FIGURE  B 
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Appendix  C 

REPORT  OF  INVENTIONS 


The  purpose  of  the  work  performed  under  this  contract 
and  reported  herein  was  the  verification  of  the  capability 
of  an  existing  remote-sensing  device  to  measure  atmos- 
pheric winds  accurately.  The  methods  for  accomplishing 
wind  measurement  had  previously  been  established,  but 
had  not  been  verified  for  laser  Doppler  velocimetry. 

Because  the  purpose  of  the  test  was  the  use  of  established 
techniques  on  an  existing  device,  no  innovations,  discovery, 
improvement,  or  invention  was  made. 
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